Density and spin response function of a normal Fermi gas at unitarity 
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Using Landau theory of Fermi liquids we calculate the dynamic response of both a polarized and 
unpolarized normal Fermi gas at zero temperature in the strongly interacting regime of large scat- 
tering length. We show that at small excitation energies the m phase (density) response is enhanced 
with respect to the ideal gas prediction due to the increased compressibility. Viceversa, the out of 
phase (spin) response is quenched as a consequence of the tendency of the system to pair opposite 
spins. The long wavelength behavior of the static structure factor is explicitly calculated. The 
results are compared with the predictions in the coUisional and superfluid regimes. The emergence 
of a spin zero sound solution in the unpolarized normal phase is explicitly discussed. 

PACS numbers: 



There is now clear experimental evidence that, in the 
unitary Umit of large scattering lengths, an harmonically 
trapped polarized Fermi gas phase separates, at very low 
temperature, into a superfluid core sorrounded by a po- 
larized normal component [H, 0, 0] . When the total po- 
larization of the gas exceeds a critical value the super- 
fluid component disappears and the gas becomes normal. 
The observed phase separation and the shape of the den- 
sity profiles in large atomic samples is theoretically well 
understood employing the equation of state of the su- 
perfluid and normal phases in the local density approx- 
imation 0, Q. These studies have revealed the crucial 
role played by the interactions in the normal phase. The 
maximum local concentration (Chandrasekhar-Clogston 
limit) of the minority (spin down) versus the majority 
(spin up) component achieved at unitarity is predicted to 
he X = ni/n^ ~ 0.45 at zero temperature, in agreement 
with the most recent experiments carried out at MIT Q . 
Larger concentrations are expected to occur if one im- 
poses an adiabatic rotation to the gas which favours the 
formation of the normal phase @. Actually in the ro- 
tating case even an unpolarized Fermi gas is expected to 
give rise, at unitarity, to phase separatation. 

The availability of the normal phase of the unitary 
Fermi gas at very low temperatures opens new stimu- 
lating perspectives, due to the expected Fermi liquid be- 
havior of this strongly interacting system. Previous stud- 
ies of the dynamic behavior have mainly focused on the 
particle (spectral) response (see, for example, 7] and ref- 
erences therein), for which relevant information is exper- 
imentally available through radio-frequency transitions 
Q. A quantitative comparison between theory and ex- 
periment is however still an open issue, especially due 
to the non trivial role played by final state interactions. 
The motion of single impurities in harmonically trapped 
configurations has also been the object of theoretical in- 
vestigation fH|. 

The purpose of this letter is to study the density and 
spin density dynamic response function of the polarized 
normal phase. With respect to the motion of single im- 
purities, the excitations investigated in the present work 
involve wavelengths much smaller than the size of the 
gas and they consequently correspond to more favourable 



conditions for reaching the coUisionless regime where the 
most interesting features of the Fermi liquid behavior 
show up. Experimentally the dynamic response can be 
measured via two-photon Bragg spectroscopy, a tech- 
nique already successfully applied to Bose-Einstein con- 
densed gases (first Bragg spectroscopic measurements 
in ultracold Fermi gases have been recently carried out 
at high momentum transfer [l3|). Differently from rf 
transitions Bragg experiments do not change the inter- 
nal atomic states and consequently final state effects are 
absent. Furthermore, by a proper choice of the detuning 
of the laser beams with respect to the atomic resonance 
one can measure suitable combinations of the density and 
spin responses. Finally, using focused laser beams, one 
can in principle measure the response locally, thereby 
providing valuable information on the uniform matter 
behavior. Theoretically the density and spin responses 
of the normal phase can be calculated using Landau the- 
ory of Fermi liquids [ll|. This theory has been so far 
mainly applied to unpolarized samples like liquid He3. It 
has been also developed to study s pin transverse excita- 
tions in slightly polarized samples jl2l . [T3| . In this work 
we consider arbitrary polarized Fermi gases. 

The basic assumption of Landau's theory is that the 
system can be described in terms of a gas of long-living 
quasi-particles interacting through a mean field. This 
assumption is usually guaranteed at sufficiently low tem- 
perature for excitations close to the Fermi surface of each 
spin component. Under these assumptions collisions be- 
tween quasi-particles can be ignored. A key feature of the 
theory is the occurrence of deformations of the Fermi sur- 
face which deeply distinguish the dynamic behavior of a 
normal liquid from the one of a superfluid. In some cases 
these deformations give rise to a new type of collective 
motion, the so-called zero sound. For sake of simplicity 
in the following we will take into account only the inter- 
action terms arising from the isotropic deformations of 
the Fermi surface, i.e. from the density modulations of 
the two spin components. In particular we will ignore ef- 
fective mass effects which arepredicted to be small in the 
polarized phase at unitarity [1, [H, [l^ [l6| . According to 
the above assumption the relevant interaction terms can 
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be derived from the interaction energy functional 



Etnt = dr eint{n\,ni) 
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The single quasi-particle (qp) Hamiltonian for each a- 
spin component (cr =| , |) will then contain an interaction 

contribution easily derivable from Eq.([T|): iJ^*^-* = Hq + 
dsint I dric where Hq is the free particle Hamiltonian. The 
quasi-particle Hamiltonian will be modified during the 
motion by position and time dependent terms, associated 
with the local changes in the density distribution of the 
two spin species. These terms have to be determined with 
a self-consistent procedure. In the linear approximation 
one has 
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where the suffix indicates that the derivatives should 
be calculated at equilibrium and 5n^ = 5n„ (r, t) are 
the space and time dependent density changes with re- 
spect to equilibrium. The dynamic response function 
is calculated by adding an external field of the form 
Vext — ao-e*'-'* '"^'^*'' to the Hamiltonian ^ and the linear 
response x<t,o-' is defined by the induced density fluctua- 
tions according to 
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By inserting result ^ into Eq.© we can calculate the 
response of the system using the free particle Hamilto- 
nian by simply replacing the external coupling with 
the effective value 
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in Vext ■ This procedure yields a self-consistency relation- 
ship for the dynamic response functions and provides the 
relevant expressions for Xo-.o-' in terms of the free particle 
responses t ^nd i (the crossed free response x^ i 
identically vanishes because there are no correlations be- 
tween spin- up and spin-down particles in the ideal gas): 



XT,t('?,'^) =Xt,t(9>^) 
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where the denominator D{q,Lu) is defined by 
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and, for simplicity, we have omitted the suffix in the 
density derivatives of Cint- The effects of the interac- 
tion is particularly crucial in the crossed xt.x response. 
As we will explicitly discuss later the sign of the rele- 
vant interaction term d'^eint/ dn^dni is negative, reflect- 
ing the tendency of the system to pair opposite spins. 
The poles of the response functions (zeros of D{q,Lu)) 
correspond to the undamped, discretized oscillations of 
the system (zero sound). They occur only in the pres- 
ence of the interaction. Useful combinations of ^ are 
provided by the symmetric (s) and antisymmetric (a) re- 
sponses Xs(a) = XT,T + Xi,i ± 2xT,i> also called density 
and spin responses. 

The concept of quasi-particle, and consequently result 
^ for the response functions, applies to small wave vec- 
tors {q <^ qpcr where qpa is the Fermi wavevector of the 
cr-specis) . In this limit the free response x" , at zero tem- 
perature, reduces to the simple form: 
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with Act = uj/qvp rr a nd Vp„ = hqp^/m. The function 
5(A) is defined by [lH 



5(A) =1-^1^1^-4^ 0<A<1 
5(A) A>1. 
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From the knowledge of the response function one can 
calculate the dynamic structure factor in each spin chan- 
nel. At zero temperature the following relation holds for 

LO>Q 
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while Scr{q,uj) = for w < 0. The function Sa-{q,c^) 
consists, in general, of a discretized peak (zero sound) 
and of a continuum of quasi-particle excitations. In the 
absence of interactions one finds S^{q,uj) — \arn/2Tr^ for 
< Act < 1 and elsewhere. 

The dynamic structure factor obeys important sum 
rules [ll[ ■ The most famous one is the /-sum rule 
/ dujujSa-{q,to) = n^q^ I2m which is exactly satisfied by 
the mean field results (O for Xct.ct, being directly related 
to the behavior of the response function at large 

w. Another important sum rule is provided by the non- 
energy weighted moment, yielding the static structure 
factor / dujScr{q,ij) — Sa{q)- This quantity is directly re- 
lated to the Fourier tansform of the two-body correlation 
function. Landau's theory provides accurate information 
on the long wavelength (small q) behavior where S^iq) 
is linear in q. 

The simplest case is the symmetric configuration (n^ — 
ni). At very low temperature the corresponding ground 
state configuration is expected to be superfluid, the 
normal phase being available only at higher tempera- 
tures where Landau's theory is not applicable. How- 
ever, as already anticipated in the introduction, if one 
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switches on adiabatically the rotation of the confin- 
ing trap, also in the unpolarized case the gas is pre- 
dicted to phase separate at zero temperature into a 
superfluid core sorrounded by a normal shell [g^, so 
that the study of the dynamic behavior of the unpo- 
larized normal gas might be relevant for future exper- 
iments. In this case the interaction coefficients deter- 
mining the response function are conveniently written as 
d'^e.ntldn^dn^ = d'^e^nt/ dn^dni = 2/i{tF/n){F^ + Fg) 

and d^e^nt/dn^dni = 2/i{eF /n){F^ - Fg) where F^''"^ 
are the Landau parameters in terms of which the density 
(s) and spin (a) response functions take the familiar form 
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The above Landau's parameters are directly related to 
the compressibility and to the magnetic suscepti- 

bility xm of the gas according to 
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where Cq = Vp/V^ and xmq are, respectively, the sound 
velocity and the magnetic susceptibility of the ideal gas. 
Since the gas is strongly interacting the evaluation of Fq 
and Fq requires a non perturbative approach. Monte 
Carlo calculations of the equation of state and of the 
magnetic susceptibility in the normal phase yield the 
values F^ ^ -0.44 ^ and F^ ^ 2 ^ at unitarity 
[20| . The behavior of the Landau's parameters of the 
unitary Fermi gas deeply differs from the one of liquid 
"^He where Fq is negative reflecting the tendency of the 
system to behave like a ferromagnet. In our case Fq is 
instead positive, reflecting the tendency of the system to 
pair opposite spins. An interesting consequence is the 
appearence of a discretized pole in the spin response, as- 
sociated with the propagation of spin zero sound, while 
the density response is Landau damped. In Fig[l] we 
show the predicted values of the density and spin dy- 
namic structure factor in the unpolarized case, using the 

2.0. The spin zero sound so- 
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-0.45 and Fq" 



lution takes place at the value uj — l.Wqvp. The different 
behavior in the density and spin channels is evident and 
shows up also in the low q behavior of the static structure 
factors: S,{q) = 1.26Soiq) and Sa{q) = OMSoiq) where 
Soiq) = i/'inq/qp is the ideal Fermi gas prediction and 
n is the gas density. The strength carried by the zero 
sound is ^ 80 percent of the total spin strength. 

In the polarized case a flrst useful description of the 
interaction effects is provided by the low x expansion 
of the equation of state. At unitarity the expansion 
takes the simple form Cmt — ~(3/5)Aei?|n| with the 
value A ~ 1 predicted by both Monte Carlo [1 HI, 
and diagrammatic calculations and where ep^ is the 
Fermi energy of the spin- up component. The expan- 
sion has been shown to provide an accurate descrip- 
tion of the equation of state up to the largest values of 




FIG. 1: (Color online) Density (dashed, black) and spin (full 
line, red) dynamic structure factor for an unpolarized Fermi 
gas in the normal phase at unitarity and T = as a function 
of A — Lojqvp. The red arrow indicates the position of the spin 
zero sound solution. The ideal gas prediction is also shown 
(dotted, blue). 



concentration achievable in the absence of rotation. In 
the small x regime the relevant interaction parameter is 

In Fig[2] we show the predictions for the density and 
spin dynamic structure factors calculated for a polar- 
ized Fermi gas at a; = 0.44. The comparison with the 
ideal gas prediction (dotted line) explicitly reveals the 
role of the interactions. Similarly to the unpolarized 
case interactions have opposite effects in the two chan- 
nels. For the static structure factors, proportional to the 
integral of the dynamic structure factor, we find the re- 
sults S.is) = 1.285o(g) and Sa{q) = 0.83S'o(g) where 
50(9) = 3/4(H-a;^/^)n|q/(7_F| = X.V&n^qjqp^ is the ideal 
gas value. The cusps in the dynamic structure factors at 
A| = 1 and A| = x^l'^ reflect the existence of two Fermi 
surfaces in the polarized case. In addition to the con- 
tinuous structure the response in the spin channel gives 
also rise to a discretized contribution (spin zero sound). 
This pole is however located too close to the continuum 
threshold A| = 1 to be physically relevant, differently 
from what happens in the unpolarized case. 

The above results for the response function hold in 
the collisionless regime of the normal phase. When colli- 
sions become important the dynamic behavior changes 
in a drastic way, being characterized by an ordinary 
(first) sound mode, giving rise to a sharp peak in 
the density response, and by a diffusive spin excita- 
tion. The first sound velocity is determined by the 
compressibility of the gas according to the thermody- 
namic relation mc^ = ndpL/dn. For the unpolarized 
case we find c = vp{l + Fq)^/^/\/3, corresponding to 
A ~ Lojqvp — 0.43, while for the polarized gas one 

finds c = {[l + x^/"^ - Ax)/{l+x)f''^ /VZ [2l|. At 
x — 0.44 the reduction of the sound velocity with re- 
spect to the ideal gas {A — 0) is still significant and 
comparable to the reduction calculated in the unpolar- 
ized case. It is also useful to compare the above pre- 
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FIG. 2: (Color online) Density (dashed, black) and spin (full 
line, red) dynamic structure factor for a normal polarized (x — 
n|/n| = 0.44) Fermi gas at unitarity and T = as a function 
of A — uj/qvpi- The ideal gas prediction is also shown (dotted, 
blue). 

dictions with the superfluid behavior of the upolarized 
gas where sound propagates hydrodynamically at the ve- 
locity c = vpil + /3)^/V\/3 with (3 ~= -0.58. In the 
superfluid the spin excitations are instead gapped with 
the threshold given by 2A where A is the single particle 
gap, of the order of the Fermi energy at unitarity. 

The transition from the coUisionless to the coUisional 
regime takes place when ujt ^ 1 where t is a typical col- 
lisional time. Using dimensional arguments one expects, 
at unitarity, the dependence H/t cx epiksT / cfY , where 



the factor originates from the PauU principle. The 
value of the coefficient of proportionality is estimated to 
be ~ 4 both in the unpolarized case [l^j and in the highly 
polarized limit [23] (in the latter case e^? is the Fermi en- 
ergy of the majority component). The conditions for ap- 
plying Landau's theory in the coUisionless regime are not 
consequently too severe, requiring frequencies satisfying 
the condition: 




The condition fuv ^ ep {ksT / ep) is in particular much 
less severe than in the case of the spin dipole oscillation 
of the trapped gas which takes place at frequencies of the 
order of the harmonic oscillator frequency (Hj- 

In conclusion we have shown that at unitarity the re- 
sponse function of a Fermi gas in its normal phase is siz- 
ably affected by the interactions and behaves quite differ- 
ently in the density and spin channels. Both mean field 
and coUisional effects are predicted to take place in ranges 
of temperatures and frequencies of reasonably easy access 
in two photon Bragg spectroscopy experiments. 
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